subroutine comput_flux1
use var
implicit none

    integer m,j,i
    real(8) f(0:gridsx,0:gridsy,4)
    real(8) local_max
    real(8) f_m,f_p
    
    call comput_Qx_m
    
    call comput_Qx_p
    
    !write(*,*)flux1(:,1,1)
    
    call L_F_local_x(Qx_m,Qx_p,f)
    do m=1,4
        do i=0,gridsx
            do j=1,gridsy
                flux1(i,j,m)= f(i,j,m)
                !if( Tstep==1 .AND. KRK==2 .AND.m==1) then
                !    write(*,*)i,j,flux1(i,j,1)
                !endif
            enddo                
        enddo  
    enddo
    
    

end subroutine